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Abstract. A general model for the cosmological evolution of the fine structure 
constant a driven by a typical Quintessence scenario is presented. We consider a 
coupling between the Quintessence scalar <f> and the electromagnetic kinetic term 
F^ U F^ V , given by a general function Bp(<fr). We study the dependence of the 
cosmological Aa(t) upon the functional form of Bp(<j)) and discuss the constraints 
imposed by the data. We find that different cosmological histories for Aa(t) are 
possible within the avaliable constraints. We also find that Quasar absorption spectra 
evidence for a time variation of a, if confirmed, is not incompatible with Oklo and 
meteorites limits. 
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1. Introduction 

Over the last few years there has been increasing interest in the possibility of varying 
the fundamental constants over cosmological time-scales. This has a twofold motivation. 
On one side, several observations point towards the existence of a smooth dark energy 
component in the universe, which could be modeled via a dynamical scalar field 
(Quintessence) pQ. In general, we expect such a cosmological scalar to couple with some, 
if not all, the terms in the matter-radiation Lagrangian, thus inducing a time variation 
of physical masses and couplings. On the other side, recent improved measurements 
on possible variations of the fundamental constants are opening up the possibility 
of testing the theoretical models to a good degree of precision over a wide range of 
cosmological epochs. It should also be mentioned that, although controversial, some 
evidence of time variation of the fine structure constant a in Quasar absorption spectra 
was recently reported [3 E]. The cosmological variation of fundamental constants 
induced by couplings with the Quintessence scalar is then worth studying in order to 
see if such a field could be responsible for a measurable effect. 

Among all the possibilities, the time- variation of the fine-structure constant is the 
simplest to study both from the theoretical and experimental points of view. In this 
paper we will then restrict ourselves to this issue. The theoretical study of a time- varying 
fine structure constant dates back to 1982 when Beckenstein j3] first considered the 
possibility of introducing a linear coupling between a scalar field and the electromagnetic 
field. More recently the Beckenstein model has been revived, generalized and confronted 
with updated experimental limits [3 El III IE] ■ The concrete case of the Quintessence 
scalar has been considered too [3 EH HH E2] ■ However, as we will see, most authors 
restrict their study to the simplest case of a linear or quadratic coupling. The possibility 
of reconstructing the dark energy equation of state from a measure of a-variation has 
also been proposed in the literature 

In this paper we will discuss a general model for the variation of the fine structure 
constant a driven by a typical Quintessence scenario. After briefly reviewing the most 
recent observational and experimental constraints on the variation of a, we will go on to 
construct the theoretical framework. In particular we will consider the case of a general 
coupling Bf(4>) (see Eq. f|2l?|) below) which includes several classes of possible functions. 
In this way we will be able to study the dependence of the cosmological variation of 
alpha, Aa(t), upon the functional form of Bp(4>) and discuss the constraints imposed 
by present data. We will find that different cosmological histories for Aa(t) are possible 
within the available constraints. We will also find that, Webb et al. data j2j 0], if 
confirmed, are not incompatible with Oklo and meteorites constraints [TH| l2Ti] . 

2. Overview of the constraints 

Comprehensive reviews about the theoretical and experimental issues connected to the 
time variation of fundamental constants can be found in Refs. jH] , ^Hj and ^E] . In the 
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following we summarize the available constraints on the time variation of a, expressed 
as functions of the redshift z (see also Fig. [IJ): 

Aa{z) _ a(z) — a , . 

a «o 
where a = a(0) is the value measured today. 

(1) The most ancient data come from Big Bang Nuclesynthesis (BBN) and give 

pa iig : 



Act 



a 



£ 10" 



10 10 - 10 s 



(2) 



(2) More recently we have the limit coming from the power spectrum of anisotropics 
in the Cosmic Microwave Background (CMB) [T9"] : 
Aa 



a 



< 10" 



(3) 



(3) From absorption spectra of distant Quasars there are more controversial data. 
Webb and Murphy's groups combined data [21 El report a 4 a evidence for a variation: 
Aa/a = (— 0.543±0.116)-10~ 5 on a cosmological time span between z = 0.2 and z = 3.7. 
This result has not been confirmed by other groups. For example, Chand et al. ^7] give: 
Aa/a = (-0.06 ± 0.06) • 10~ 5 for z = 2.3 - 0.4 and Aa/a = (0.15 ± 0.43) • 10~ 5 for 
z = 2.92 — 1.59. We have chosen to be conservative and will consider a limit based on 
the last two results, which is consistent with zero variation: 
Aa 



a 



£ 1(T 



3-0.4 



(4) 



(4) From the analysis of the ratio Re/Os in meteorites dating around 4.56 billion 
years ago it is possible to compute 187 Re half-life, which gives [TH] : 
Aa 



a 



£ 10' 



0.45 



(5) 



(5) From the Oklo natural nuclear reactor that operated 2 billion years ago in 
Gabon, we also have [T51 l2Tij : 
Aa 



a 



£ 10 



-7 



0.14 



(6) 



(6) We then have limits coming from laboratory measurments which constrain the 
present rate of change of a. Comparing atomic clocks, which use different transitions 
and atoms, what is obtained is l2"Tl: 



< 10 



yr- 1 







(7) 



where the dot represents differentiation w.r.t. cosmic time. 

(7) In addition to the limits discussed above, there is a constraint coming from 
indirect violation of the Weak Equivalence Principle (WEP), parametrized by the Eotvos 
ratio 

\ a l — 02 1 



V 



ai + a 2 \ 



(8) 
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Log of alpha variation 
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Figure 1. The experimental constraints (l)-(6) discussed above are summarized in 
the picture: log|Aa/a| is plotted as a function of the redshift z. On the right-hand 
side we zoom on z ^ 10. The grey areas are those excluded by present data. 



where a\ e 02 are the accelerations of two different test bodies in the Earth gravitational 
field. These constraints come from the fact that nucleon masses get electromagnetic 
corrections from quark-quark interactions. As extensively discussed in j22j, the leading 
term of the electromagnetic contribution comes from the electrostatic energy of the 
quark distribution, which is proportional to a. The corrected masses, to leading order 
in a, can then be written as [TT| I22j: 

m p = m + a B p 

m n = m + a B n (9) 

where p, n stand for proton and neutron and B p = 0.63MeV/a ? B n = — 0.13MeV/a - 
If we suppose that a = a (</>), then we will induce a 0-dependence on the nucleon masses: 

5m n = B n Sa ; Sm p = B p 5a . (10) 



If we define 



we get an indirect violation of the equivalence principle induced by the 'fifth-force' 
mediated by the scalar field 



where: 



M 2 

V - ^2 (Rn9n + Rp9 P ) (&Rn9n + ARpgp) (12) 

Rf = J 1 ' = ~ 0.5 AR t ee |nM ~ ni ' 21 ~ 0.06 - 0.1 , (13) 
+ nZ n n + n p 

and fh ~ 931 MeV is the atomic mass unit. The suffix E refers to the Earth, while 
1 and 2 refer to two test bodies having equal mass but different composition. From 
Eqs. (JTIJ)-(JT2J) we see that any model of a- variation will induce a characteristic g Ptn 7^ 
and hence WEP violation: while the first two factors in Eq. ()12j) are universal and 
depend on the Earth composition, the third term is not zero if and only if g p n 7^ and 
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the test bodies have different composition in neutrons and protons. The current limits 
on WEP violations impose 



77 < 10 



-13 



(14) 



3. The theoretical framework 

Following Olive et al. the most generic action involving a scalar field, the Standard 
Model fields and an hypothetical Dark Matter particle Xi can be written as 



S 



16ttG 
- / d^x 



+ / d A x 



+ 



f d 4 x^g R + J d 4 

—g B F {<P)F, V F^ - 1 J d'x^g B^F^F^ 

j 

d^xyf^ [x$X ~ B x ((f))m x x T x] (15) 

where p= j^D^ and for the electromagnetic term, for example, = <9 M — ie^A^. The 
index i — 1, 2, 3 refers to the SU(3) gauge group of the Standard Model and j runs over 
the various matter fields. 

The form of the action (JT5|l follows from supplying ^-dependent factors to all mass 
and kinetic terms to the standard Lagrangian (which would have all Bi — 1). In general 
we would expect that all of the l?j(0)'s are switched on, if not forbidden by any symmetry 
principle. However, the theoretical treatment of the full Lagrangian is very cumbersome 
and so the coupling functions Bi(<f>) are usually switched on one at a time. In this way 
one can also disentangle the effects due to each single term. In this paper we want to 
focus on the fine-structure constant a, and so we will keep only Bf(4>) 7^ 1 and set all 
the other functions equal to 1. 

The relevant part of the action for the effect we are going to study is then 



S 



16vrG 
d 4 X\ 



d 4 X\f^g R + J d A x\ 
~g B F (4>)F^F^ 



-g 



-dyo^ - v(<f>) 



which allows to define an "effective" fine structure constant 

a 



a(t) 



(16) 



(17) 



BfW)) 

where is the value measured today. From (|T7|) we obtain the relative variation relevant 
for each cosmological epoch 

Aa _ a(t) - a _ 1 - B F {<t>{t)) 

a ~ a ~ B F ((f)(t)) ^ ' 

It can immediately be seen that, depending on the cosmological evolution of <j)(t) and 
on the functional form of B F ((p), the fine structure constant a could in principle have 
had many possible histories during the life-time of the universe. What possibilities are 



Cosmological evolution of Alpha driven by a general coupling with Quintessence 



6 



allowed by a general coupling Bp(<fi) within the available observational constraints is 
then worth studying. 

The relevant equations governing the cosmological evolution in a flat universe are 
the following 

! = -^£(l + **>A (19) 

ff 2 = (-)=3^I> PO) 

• dV 

+ 3F0+ — = 21 
d(p 

where i = m, r, <fr runs over the matter (including dark matter), radiation and scalar 
components. The relevant equations of state are w m = for matter, w r = 1/3 for 
radiation and as defined in Eq. ()22|). It is important to note that the evolution 
equation of the Quintessence scalar (}2~T]) does not depend on Bp or its derivatives. This 
is due to the fact that the statistical average of the term F^F^ over a current state of 
the universe is zero. So the only term that drives <fi during the cosmological evolution 
is the potential V(4>). 

Since we are working under the hypothesis that the scalar field <fi in Eq. (|2*T|) is 
the Quintessence scalar, we should also impose the additional constraints coming from 
Quintessence phenomenology. In particular we choose a runaway potential which goes 
to zero as far as the field <ft rolls to infinity, in accordance with the observational data. 
It is also required that the scalar dynamics gives the correct value for the equation of 
state 

(S -°' 7today) (22) 

The most general form for the Quintessence potential involves a combination of a power- 
law and exponential terms [21] . For the purpose of this paper, however, we will consider 
the simplest case of an inverse power-law potential V(<f>) = M 4+n 0~ n , which gives a 
late-time attractor equation of state = — 2/(n + 2) during matter domination [25]. 
The potential should also be normalized in order to give the correct energy density today 
(p9 ~ V((f>) ~ 2/3 p°): this sets the mass scale M. In what follows we will choose n = 1 
in the potential in order to have the correct attractor equation of state, and so obtain 



M ~ ^2/3 p° c M p . 

We have checked that choosing different Quintessence potentials gives a 
subdominant effect on the cosmological variation of a, with respect to changing the 
coupling function Bf{4>)- In what follows we will then fix V(<f)) = M 5 /(p and study the 
effect of different Bp(<f>ys. An interesting study, which is complementary to what is 
done here, is that of Ref. [H] where the effect induced by different Quintessence models 
on the cosmological Aa is examined in detail, while keeping the function Bp{ip) fixed. 

In order to be as general as possible we will consider a function Bp{(j>) which is a 
combination of different possible behaviors and characterized by a set of four parameters 
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that are allowed to vary freely: 

Bf{4) = (J^J [1 - C(0 - <Po) q ] . (23) 

This choice is not motivated by a specific theoretical model, but is rather a working 
tool for obtaining different functional forms of Bp{<f) and thus cosmological histories 
of a, according to Eq. (fTSj). We have chosen a combination of functions (power-law, 
polynomial, exponential) that can be switched on and off at will (depending on the 
values of the parameters e, £, r and q), thus giving rise to a variety of possibile .Bf(0)'s. 
In this way we can carry on a unified discussion of a number of different models of a 
variation. 




Figure 2. Evolution of the energy densities (left) and scalar equation of state (right) 
for a quintessence model with potential V — l/4> and initial conditions p™/Pc = 10 30 
at z = 10 10 . The dot-dashed line represents the energy density of radiation, the 
dotted line the energy density of matter, the green dashed line the energy density 
of quintessence and the red solid line the attractor. All of the energy densities are 
expressed in units of the present critical energy density p° c . 



4. Cosmic evolution of a 

We have numerically solved the cosmological equations (fE?J) -fl2T ]) and then plotted the 
resulting cosmological history of Act for various classes of functional forms of Bf{4>), 
according to Eq. (jl8j) . As already mentioned, for illustrative purposes we have chosen a 
scalar potential V — 1/0 and initial conditions j ' p° c = 10 30 at z = 10 10 . Fig. |2 shows 
the corresponding evolution of the energy densities and of the scalar equation of state 
parameter. 

Linear coupling 



The simplest case is given by the choice e = r = and q = 
Eq. fl21 : 

£fc(0) = l-C(0-0o) 



1 for the parameters in 
(24) 
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This case corresponds to the original Beckenstein proposal jlj, which however did 
not supply a potential to the scalar field§. Copeland et al. jH] give a comprehensive 
discussion on various Quintessence models linearly coupled to the electromagnetic field, 
but assuming Webb et al. data [21 IH] to be correct and imposing on Aa(t) to agree with 
that measure. 

In our case, the resulting Aa, as defined in Eq. (JTSj) . is plotted in Fig. El We tried 
a number of different values for (, in order to verify in which cases all the available 
experimental constraints were simultaneously satisfied. We found that they are all 
respected for ( < 0.6 • 10~ 6 . With this choice, the constraints on the violation of 
equivalence principle and the constraints derived from atomic clocks are automatically 
satisfied: 

- ^—9i . . , ^ — i a & 

V 



4 ■ 1(T 21 < 10 -13 



4 • 10" 17 < 10 



-15 



(25) 



Log of alpha variation 




-Log z+1 



Figure 3. The logarithm of |Aa/a| is plotted as a function of Log(z + 1) for 
B F (4>) = 1 - C(<£ - 4>o) with C = 0.6 ■ 1CT 5 (dotted line), C = 0.6 ■ 10~ 6 (solid line) and 
C = 0.6 • 10 -7 (dashed line). On the right-hand side we zoom on z ^ 10. Only the 
curves not overlapping the grey areas are phenomenologically viable. 



Polynomial coupling 

A slightly more complicated case is given by the choice e = r = 0, allowing the exponent 
q to be > 1: 

5 F (0) = l-C(0-0o) 9 (26) 
The case of a quadratic coupling (q = 2) was considered in Ref. [TU], but with the 
additional assumption of a proportionality relation between B F (<f)) and V (</>). 

We have found that the data do not impose any upper limit on the exponent q 
and that increasing q makes it possible to reduce the fine-tuning in (. For example, 
choosing ( = 10~ 4 the experimental limits are respected for q = 6 and the constraints 

§ To be precise, Beckenstein actually invoked an exponential coupling, which however is practically 
equivalent to eq. (|24Jl due to the smallness of 
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on the violation of equivalence principle and the constraints derived from atomic clocks 
satisfied by many orders of magnitude. In Fig. 0]we plot Log |Aa/a| for £ = 10~ 4 with 
q = 3, 6 and 9, as function of red-shift. 



Log of alpha variation 




-Log z+1 -Log z+1 

Figure 4. The logarithm of |Aa/a| is plotted as a function of Log(z +1) for 
Bp{4>) = 1 — C{4> — (f>o) q with £ = 10~ 4 and q = 3 (dotted line), q = 6 (solid line) and 
q — 9 (dashed line). On the right-hand side we zoom on z <J 10. Only the curves not 
overlapping the grey areas are phenomenologically viable. 



As already mentioned, increasing the exponent q we can do even better. For 
example, with q = 17 the experimental constraints are satisfied even for ( = 1, as 
illustrated in Fig. El It should be emphasized that among all the possibilities considered 
in this paper, this choice of the parameters appears to be the most natural of all. A 
notable feature is also the fact that the value of Aa is enhanced in the past, with 
respect to the other cases, becoming closer to the observational limits, while falling off 
very steeply in recent times. 



Log of alpha variation 
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Figure 5. The logarithm of |Aa/a| is plotted as a function of Log(z + 1) for 
Bp((f>) = 1 — C(0 ~ 4>o) 9 with ( = 1 and q = 17 (solid line). On the right-hand 
side we zoom on z ^ 10. Note that all the experimental limits are satisfied without 
any fine-tuning in the parameters of the function Bp{4>)- 
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Power-law coupling 
With the choice ( = 

B F {4>) 



we obtain the following coupling function: 




(27) 



In this case, it is necessary to fine-tune the exponent e in order to satisfy the data, 
due to the smallness of <fi in the early universe. Keeping e of order one would violate 
even the constraints from BBN. We found that the experimental limits are respected for 
|e| < 4- 10 -7 . In Fig.lHlwe plot Log |Aa/a| as a function of red-shift for different choices 
of e. Note that the sign of Aa depends on the sign of e. With the choice e = 4- 10~ 7 , the 
constraints on the violation of equivalence principle and the constraints derived from 
atomic clocks are automatically satisfied: 

' a 



rj ~ 4- 10~ 21 < 10 



,-13 



Q'o 



4 • 10~ 17 < 10~ 15 yr" 1 



(2f 



Such a small exponent might look quite unnatural, however Eq. (|27|) for e 1 is 
equivalent to: 




B F {<j>) = 1 + eln 



In this way the fine tuning is moved from the exponent to the coefficient. 



(29) 



Log of alpha variation 
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Figure 6. The logarithm of |Aa/a| is plotted as a function of Log(z + 1) for 
B F {4>) = (j^y with e = 4 - 10" 6 (dotted line), e = 4- 1CT 7 (solid line) and e = 4 - 10" 8 
(dashed line). On the right-hand side we zoom on z ^ 10. Only the curves not 
overlapping the grey areas are phenomcnologically viable. 



Exponential coupling 

The choice e = ( = of the parameters in 

B F (<P) = e - T ^- M . 

In this case, if r ^ 1 it is not possible to 
Depending on the sign of r, the resulting | 



(P3jl gives: 

(30) 

satisfy all the constraints at the same time. 
Aa/a\ becomes too large in the early or late 
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universe. For r« 1, instead, the coupling function becomes equivalent to the linear 

case: B F {<j>) = e^"^ ~ 1 - r(0 - <f) ). 



Linear and power-law coupling combined 

Now let's consider two factors in (|23|) with q = 1 and r = 0. 

Bf{(P)= (£) ( i_c ^~^ • (3i) 

For an arbitrary choice of ( and e, the resulting Aa is similar to the linear coupling or 
power-law coupling case, depending on which factor dominates. It is instead interesting 
to consider the case ( = 7e in which the two factors are of the same order of magnitude. 
If 7 > 0, the two factors can contribute in an opposite way and it is easy to obtain 
Aa ~ also at some time in the past. For example, with the choice e = 2.4 ■ 10~ 6 and 
7 = 2.2 we obtained the behavior plotted in Fig. [?l in which we have varied 7 of 10%. 
For the choice e = 2.4 • 10~ 6 , 7 = 2.2 the constraints on the violation of equivalence 
principle and the constraints derived from atomic clocks are automatically satisfied: 



V 



2 • 1(T 20 < 10" 13 



a 



Q'o 



9 



1(T 17 < 1(T 15 yr" 1 



(32) 



Log of alpha variation 
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Figure 7. The logarithm of |Aa/a| is plotted as a function of Log(z + 1) for 
B F {4>) = (^) £ (1 - 7e(0 - M) with e = 2.4 • 1(T 6 , 7 = 2.2 (solid line) and 
7 = 2.2 ± 10% (dashed and dotted respectively). On the right-hand side we zoom 
on z ^ 10. Only the curves not overlapping the grey areas are phenomenologically 
viable. 



Power-law and exponential coupling combined 

Since, as already discussed, the exponential coupling function case is equivalent to the 
linear one, this possibility falls within the previous example. 
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Polynomial and exponential combined 

Now let's consider e = and g ^ 1. The case q = 1 is not interesting since the two 
factors would be almost equivalent and the behavior corresponding to the linear case. 
Let's choose then, for example, q = 6: 

S F (0) = (l-C(0-0o) 6 )e-^ o) (33) 

For recent times (z < 1) the exponential coupling e~ T ^~^ ~ 1 — r((f) — 4>q) dominates, 
while in the past the two terms can be of the same order and, due to q being even, cancel 
at some time. This is shown in Fig. |H1 The constraints on the violation of equivalence 
principle and the constraints derived from atomic clocks are satisfied. 



Log of alpha variation 
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Figure 8. The logarithm of |Aa/a| is plotted as a function of Log(z + 1) for 
B F ((j>) = (1 - C {4> - 0o) 6 ) e-^-^ with r = 0.6 • 1CT 6 and ( = 2 • 1(T 4 (dotted line), 
C = 3.2 • 10~ 5 (solid line) and ( = 5 • 10~ 6 (dashed line). On the right-hand side we 
zoom on z ^ 10. Only the curves not overlapping the grey areas are phenomenologically 
viable. 



5. Summary and Conclusions 

In this paper we have carried out a comprehensive study of the cosmological variation of 
the fine-structure constant a induced by the coupling of the electromagnetic field with 
a typical Quintessence scalar. We have considered a variety of functional forms for the 
coupling function Bp ((ft), obtainable from a general expression (see Eq. (J221)) depending 
on four parameters. 

We have found that very different cosmological histories for Act are possible, 
depending on which parameters are switched on. For example, we can produce a Ace 
which is well below the observational constraints in the early universe and just within 
the experimental limits in recent times (linear coupling case). But also the converse 
is possible, if we choose a polynomial coupling. In particular, the behavior at small 
red-shift can be qualitatively very different depending on the model we choose: sharply 
decreasing in the polynomial coupling case or mildly decreasing with the power-law 
coupling. 
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By combining different functional forms, a notable feature emerged. In some cases 
it is possible that the scalar dynamics drives Aa to a zero at some time in the past, 
thus inverting the slope of its cosmological evolution. This happens for all the combined 
cases discussed here. 

It is also worth remarking that in our parameter space span we found solutions with 
extremely reduced fine-tuning and still compatible with the available constraints. This 
is the case of the polynomial coupling with exponent q > 15 which lifts the fine-tuning 
of the coefficient ( (usually constrained to be < 10 -6 ) to order 1. 

It should also be emphasized that, while in the literature the result by Webb et 
al. [21 E! is usually said to be incompatible with the Oklo limit 20|j, we found that this 
is not always the case. For example, in the polynomial coupling case it is possible to 
obtain several examples with a Aa which matches the Quasars data and at the same 
time respects the Oklo bound. 
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